Using quasars as standard clocks for measuring cosmological redshift 
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We report hitherto unnoticed patterns in quasar light curves. We characterize segments of quasars' 
light curves with the slopes of the straight lines fit through them. These slopes appear to be directly 
related to the quasars' redshifts. Alternatively, using only global shifts in time and flux, we are able 
to find significant overlaps between the light curves of different pairs of quasars by fitting the ratio of 
their redshifts. We are then able to reliably determine the redshift of one quasar from another. This 
implies that one can use quasars as standard clocks, as we explicitly demonstrate by constructing 
two independent methods of finding the redshift of a quasar from its light curve. 

PACS numbers: 



Introduction. To probe the largest distances in our 
universe, we seek the brightest sources. Quasars are cer- 
tainly among these, and they are numerous. Unfortu- 
nately quasars exhibit large dispersion in luminosities at 
all wavelengths. This makes them unusable as standard 
candles for measuring cosmological distances. In an 
interesting correlation between emission line equivalent 
width (the ratio of integrated line flux over local con- 
tinuum flux density) and rest-frame ultraviolet luminos- 
ity was observed. Namely, C IV 1549 A emission-line 
equivalent width in quasars decreases with increasing UV 
continuum (1450 A) luminosity. Since flux ratios are 
distance-independent, this anti-correlation could allow 
for quasars to be used as standard candles. However, the 
large dispersion in this anti-correlation gives very poor 
distance calibrations compared to other standard can- 
dles. It is therefore worthwhile to look for other correla- 
tions that might make quasars viable distance standards. 

Identifying the patterns in quasars' light 
curves. In this paper, we analyze the V- 

band light curve data from MACHO spectroscopically- 
confirmed quasars behind the Magellanic Clouds [2|- 
These light curves and spectra are available at 
http : 1 1 www. astro. yale.edu I rngeha I MACHO j\ Most 
analyses of quasar light curves have considered the vari- 
ation of their absolute magnitude 0-0], whereas we will 
use the flux, which is more physical. The relation is 

m = -2.51og 10 (F/F°). (1) 

Here F is the flux, F° is a V-band-spccific normalization 
constant (3.636 x 10~ 20 'er gs /cm 2 /sec/ hz), and m is the 
quasar's absolute V-band magnitude. 

To study intrinsic patterns of quasar activity it is best 
to transform to the quasar's rest frame, i.e. at minimum 
we need to rescale the observed time by (1 + z)^ 1 . In 
Fig. [TJ we plot the flux, F, as a function of quasar-rest- 
frame time for selected quasars from among the MACHO 
sample. Two distinct patterns of time variation are ap- 



parent as linear trends in these rest-frame light curves. 
To emphasize these trends we mark them with parallel 
lines of two distinct slopes. Shorter time-scale trends 
(marked by the steeper thin black lines) are superim- 
posed on the longer time-scale trends (thick green lines) . 
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FIG. 1: The sample of quasars with different redshifts. Time 
is transformed to the quasar rest frame (i.e. divided by (z + 
1)). A pattern emerges, some parallel lines (same slopes) 
appear. We mark them with the thin black lines. A possible 
shallower, long-term trend is indicated by thick green lines. 

We focus on data with consistent observations over an 
interval longer than 90 days, since that allows for more 
reliable identification (or rejection) of the observed pat- 
terns. We pick several straight line light curve segments 
from 13 quasars and shift each segment so that it is cen- 
tered at the origin (t = 0,F = 0). If the slope of the 
segment is negative, we multiply it by — 1, since we are 
interested only in the absolute value of the slope. We 
then fit a straight line through the data collected from 
the corresponding segments of all 13 quasars and find 
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FIG. 2: The blue points are data points chosen from 13 
quasars: 6.7059.207; 13.5962.237; 9.4882.332; 11.8988.1350 
208.16034.100; 207.16316.446; 68.10972.36; 75.13376.66 
5.4892.1971; 59.6398.185; 13.5717.178; 13.6805.324: 
1.4537.1642. The time and flux are shifted so that 
each segment is centered at the origin (t — 0,F — 0). The 
red line is the least-squares fit for the data. The slope is 
k = 6.47 ± 0.07 x 10~ n units/day. 



the slope, (see Fig. We assume that all data points 
have the same weight in order to avoid a single segment 
giving a dominant contribution. The slope k from the 
least-squares fit is 



k = 6.47 ± 0.07 x 10~ n units/day. 



(2) 



where the "unit" is F°. If this pattern of parallel slopes in 
quasar light curves is not a coincidence, it would indicate 
that one can use quasars as distance standards, e.g. like 
Type la supcrnovae. 

Finding the redshift of a quasar from its light 
curve - Linear fit. We now illustrate the procedure of 
determining the redshift of an unknown quasar from its 
light curve. We choose a quasar that was not included in 
our sample of 13 quasars. We identify five segments from 
its light curve that appear straight and parallel to one 
another, as shown in Fig. [3] (the number of segments is 
not crucial - the more the better) . We take each of those 
segments, discarding the rest of the light curve, and shift 
them so that each segment is centered at (t = 0,F = 0). 
Note that the time axis is in observer time, not quasar 
rest-frame timesince the redshift is "unknown". We fit 
the data collected from all the segments (all belonging 
to the same quasar) with a single straight line and find 
the slope, as shown in Fig. @] In this case, the slope is 
1.94±0.16 x 10~ n units/day. Comparing this slope with 
that in Eq. (J5|) for the 13 quasars rest-frame light curves, 
we can calculate the redshift of that particular quasar 



6.47 ±0.07 x 10" 11 
1.94 ±0.16 x 10- 11 



- 1 = 2.34 ±0.27 



(3) 



of that quasar, z = 2.32. 

The same procedure can be repeated for the slower 
variations (marked by the thick green lines in Fig. [1). 
The results are again consistent within the statistical er- 
rors, though the error bars on the slope are larger - pre- 
sumably because of the "noise" from the fast variations 
as well as from other less coherent variations. 




FIG. 3: The light curve from quasar 9.5484.258 which was 
excluded from the step shown in Fig. 2. We choose for the 
linear fits those segments marked by shaded rectangles. 
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Indeed, this value is consistent with the known redshift 



FIG. 4: Graph similar to Fig. O but only one quasar is in- 
cluded whose redshift we want to calculate. The best fit slope 
is 1.94 ± 0.16 x 10 -11 units/day, which implies a redshift of 
2.34 ± 0.27. This agrees with the measured value z = 2.32. 



Therefore, if one can identify an oscillation mode that 
the particular segment of the light curve belongs to, one 
can readily calculate the redshift from the light curve. 
Note that this implies that the time dilation (which is 
simply a counterpart of the redshift) is included in the 
quasar light curves, in contrast with conclusions in [H . 

Finding the redshift by matching the light 
curves. Rather than identifying the segments through 
which we can fit a line as we did above, we could try to 
match the whole (or significant portions of the) quasar 
light curves. We do not expect to match the curves from 
an arbitrary pair of quasars as many different classes 
of quasars exist. Fig. [5] shows the observed V-band 
light curves of two quasars of nearly identical redshifts: 
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206.17052.388 (z = 2.15) and 25.3712.72 (z = 2.17). 
The similarity is not apparent at all. However, for com- 
parison, we can shift the time and flux origin of quasar 
25.3712.72 (and flip it): 

t n = t + 10 F n = 6.3 x 1(T 8 - F (4) 

t n and F n are the time and flux after the shifts and flip. 



206.17052.388 
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FIG. 5: Light curves of two quasars. The black and red lines 
belong to quasars 206.17052.388 and 25.3712.72 respectively. 
The time has been normalized to the quasar's rest frame. 
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time shift (At) which just resets the initial time and has 
no special physical meaning. The second one is a global 
offset of the flux intensity (AF). Our example in Figs. [5] 
and [5] shows that F n (t) — > — F n (t) is also a necessary 
transformation to allow. Finally, we must rescale the 
time coordinate by a factor an = -f+f^T- ^ i g °f course 
this factor a\ 2 that we are most interested in obtaining. 
Therefore, the complete set of allowed global transforma- 
tions of the flux is 



F n (t) ->• ±F n (a 12 (t - At)) - AF 



(5) 



Applying transformations in Eq. (J5J) to 5.4892.1971, we 
can find a very good match with 11.8988.1350, as shown 
in Fig [71 We also get an approximate value of ai 2 . Once 
we have matched the curves, we can identify the region 
where the light curves have the best overlap. It is that 
region that we will use to fit the data statistically. 
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FIG. 7: The light curve of the quasar 54892.1971 (z = 1.58) 
is transformed using only global transformation as in Eq. {5} 
in order to match with 11.8988.1350 (z = 0.33). Matching the 
curves with significantly different redshifts perhaps eliminates 
contamination of space as the explanation of similarities. The 
gray rectangle represents the region of the greatest overlap 
which we will use for our fit. 
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FIG. 6: The same quasars as in Fig. with the red one 
transformed according to Eq. (TjTJ). The match is now apparent. 

Fig. |6] shows that the shifted light curves overlap with 
one another remarkably well. 

We now define a procedure of finding the rcdshift of a 
quasar from its light curve. To do this, we start with the 
observed quasar light curves, for example 11.8988.1350 
(z = 0.33) and 5.4892.1971 (z = 1.58) We will keep 
the original data from 11.8988.1350, and manipulate the 
data from 5.4892.1971 to obtain a match. As in the 
previous examples, we will use four kinds of global trans- 
formations to shift the light curves. The first one is the 



The main goal is to find the value of ai 2 = l+ll wn i cn 
gives the best match to the data. From FigUJ we identify 
the region where the light curves have the best match 
(gray area). We extract the light curve data from these 
regions and transform one of them using At, AF and a. 
We then combine the data from these two light curves 
into a single light curve, Flux(At, AF, a^', ti), according 
to their new time sequence. 

To check if a given transformation gives a good match, 
we fit Flux(Ai, AF, ; £>) with a quartic polynomial: 

S(ti) = a 4 tf + d3tj + a 2 t\ + aiU + a (6) 

Fig. E] shows the best fit, for the three-parameter (At, 
AF and a\ 2 ) At- We have minimized 

X 2 ^(Flux(^)-<^)) 2 (7) 

i 
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with respect to At, AF and a\2- The standard deviation 
is a 2 = x 2 1 (N — 1), where N is the total number of data 
points. Fig.|9]shows x 2 m units of a 2 as a function of at\2- 
(We have minimized with respect to At and AF.) Within 
one cr (or 68% confidence level) a = 0.508^;^f. This is 
consistent with the actual measured value of 0.516. 
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FIG. 8: The red line is the combined light curve 
Fhix( At , AF, y+j 1 , ti) and the black line is the fitting function 
S(ti) (in our case quartic function) 
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quasars behind the Magellanic Clouds [2j, we observed 
patterns in quasar light curves that have previously gone 
unnoticed. We analyzed the light-curves in two ways. 
First, we characterized segments of the light curves by 
the slopes of straight lines through them. These slopes 
appear to be directly related to the quasars' redshifts. 
This allowed us to formulate a method for determining 
the redshift of an unknown quasar from its light curve. 
The results match the known values extremely well. This 
technique appears to allow us to obtain the redshifts of 
quasars with such linear trends within a few percent. We 
also formulated an alternative method for determining 
the redshift that does not rely on a linear fit. Matching 
the segments of two quasars light curves, we were able to 
fit for the redshift ratios of the quasars, again within a 
few percent. 

These techniques suggest that similar patterns shared 
by different quasars may allow them to be used as stan- 
dard clocks (or candles) to quantify luminosity distance. 

We performed our analysis for V-band light curves, 
though a similar procedure could be carried out for other 
wavelengths. We currently do not have a theoretical ex- 
planation of this effect. We would not want to specu- 
late much on the possible explanation since the physics 
of these objects is poorly understood. Hcuristically, if 
the frequency, /, and the corresponding amplitude, A, 
of the oscillation mode satisfy f A =constant, and one 
looks at the sine- wave oscillations, then a constant slope 
Asin(ft) ~ Aft would appear for small ft. Alter- 
natively, it could happen that particular quasi-periodic 
quasar oscillations described in Q (see also @ where 
similar objects are studied) are behind this effect. For 
related studies see also (9h11|. To identify the physics of 
the pattern we discovered clearly requires further inves- 
tigation. 

Regardless of its theoretical explanation, this observed 
effect suggests that one might be able to use quasars as 
distance standards. It will be important to extend the 
study to other quasars in order to find the dispersion of 
this effect. This will require a larger sample than the 14 
high quality quasar light curves, with continuous cover- 
age over at least 90 days, that are currently available. 



FIG. 9: x in units of °" for fitting Flux(A£, AF, i±£L, U) 
with the quartic function, x 2 nas a minimum at the redshift 
ratio of 0.508. The measured value is 0.516. 

Conclusions. By studying the data from MACHO 
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